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CALCULATION  OF  THE  SPATIAL  DIFFUSION 
OF  SMALL  PARTICLES  IN  AN  EXPONENTIAL  ATMOSPHERE 

by  Milton  M.  Klein,  Kwang  Yu,  and  Richard  H.  Harrison 

1 .  Introduction 

The  problem  of  general  diffusive  motion  of  small  particles  through 

an  atmosphere  of  variable  density  has  been  treated  recently  by  several 

investigators.  A  numerical  solution  for  the  case  of  an  exponential 

atmosphere  has  been  made  by  Banister  and  Davis'.  An  analytic  solution 
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motivated  by  their  results, has  been  obtained  by  Granzow  in  the  form  of 

an  infinite  series  of  Laguerre  polynomials.  In  one  analysis  of  the  F£ 

layer  of  the  ionosphere  as  a  charged  system  in  which  diffusion  occurs  in 

3 

the  ambipolar  limit,  Yonezawa  has  obtained  a  solution  which  may  be 
applied  to  neutral  particles  whose  mass  is  half  that  of  the  ambient 
particles.  These  analyses  have  been  limited  to  the  case  of  one-dimen¬ 
sional  motion.  Recently  Klein  and  Yu"  have  obtained  analytic  solutions 
for  both  three-dimensional  and  one-dimensional  motion  in  atmospheres 
having  exponential  and  parabolic  density  distributions.  The  numerical 
work  was  confined  to  the  one-dimensional  problem  for  the  exponential 
atmosphere  An  extension  of  the  calculation  to  the  three-dimensional 
problem  is  presented  here  for  the  exponential  atmosphere. 
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2 .  Ana  lysis 


The  diffusion  of  small  particles  in  an  atmosphere  whose  density 
varied  exponentially  with  altitude  z  was  studied  in  detail  in  Reference 
A.  The  diffusion  coefficient  D  was  taken  in  the  form 


D  z/H 

—  =r  n 

D 

o 


(1) 


where  is  a  reference  value  of  D  at  z  =  0,  and  11  is  the  scale  height. 
An  analytic  solution  for  the  three-dimensional  density  distribution 
n(r,t)  was  obtained  in  the  form1* 
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where  n  is  the  total  number  of  particles  r  -  e*Z'2H,  r  =■  e*zo/2H> 
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‘  +  (a  - 1 )  ,  and  a  is  the  ratio  of  the  mass  of  a  foreign  particle  to 

that  of  an  ambient  particle.  In  Equation  (2)  the  unit  of  length  is  H 
while  the  time  t  is  measured  in  units  of  H  /D  (dimensionless  values 


of  r  and  t  are  primed). 


Numerical  calculation  of  the  density  .’i  requires  the  evaluation  of 
the  definite  integral 
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occurring  in  Equation  (2).  Because  of  the  unusual  requirement  of  inte¬ 
gration  of  a  Bessel  function  with  respect  to  a  variable  which  depends 
upon  the  order  of  the  Bessel  function,  an  exact  determination  of  the 
integral  N  does  not  appear  possible.  An  approximate  analytic  solution 
for  the  late  times  may,  however,  be  obtained.  For  larlge  values  of  t  the 

parameter  X  «  1  and  I  (X)  may  be  expanded  in  a  powep  series 
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and,  upon  retention  of  the  first  two  terms,  Equation 
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(3)  becomes, 
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The  solution  of  Equation  (3)  therefore  depends  upon  Expressing  the  gamma 
function  (3+2)  in  a  suitable  form.  The  fact  that  f]or  large  x  the 
asymptotic  form  of  the  gamma  function  is 

r*  Jx  iosx !  *x 
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suggests  that  a  suitable  form  for.  r(p+2)  that  yields  a  tractable  integral 
which  is  reasonably  accurate  is 

’  (3+2)  *  c  3  ebl<  ,  (7) 

where  c  and  b  are  constants  which  can  be  evaluated  by  matching  Equation 
(7)  to  the  exact  gamma  function  at  two  suitable  points.  Equation  (7) 
is  most  accurate  at  large  3,  ie.,  for  large  a,  but  ^ends  to  give  too 
large  values  near  a  =■  1  Utilization  of  the  integral 
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where  a  «  a-1  yields  for  the  integral  N 


where 
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The  integral  N  therefore  exhibits  an  exponential  behaviour  e  ^  2  JJ  for 
large  o. 


For  the  case  a  «  l,  where  0  *  ,  ,  the  exponential  vanishes  and  the 
integral  may  be  represented  as  follows: 
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3 .  Results  and  Discussion 


Some  typical  results  obtained  from  this  investigation  are  presented 

in  Figures  1  through  9.  In  Figures  1  through  3  we  have  plotted  the 

relative  density  n/n  (n  *  n  at  a  *=  0)  as  a  function  of  j  for 
max  max 

several  typical  values  of  a  at  600,  400,  and  150  km.  In  all  cases  we 
note,  as  anticipated,  that  the  heavier  particles  exhibit  a  more  rapid 
drop-off  in  density,  regardless  of  the  altitude,  than  do  the  lighter 
particles.  It  should  not  be  inferred,  however,  that  the  absolute  density 
at  a  particular  location  for  a  heavy  particle  is  always  below  that  for  a 
light  one. 

Absolute  density  contours  in  the  z  -  j  plane  for  a  «  1,  2,  5  and  10 

are  represented  in  Figures  4  through  9  at  time  t'  «  10  \  We  note  that 

as  the  altitude  becomes  lower  the  densities  for  the  heavier  particles 
(a  5,  for  instance)  are  greater  than  those  for  the  light  particles 

(a  »  1,  for  example)  A  comparison  of  the  a  «  1  and  a  =  5  contours 

(Figures  3  and  7)  indicates  the  very  rapid  spreading  out  of  the  lighter 
particles  and  the  extremely  small  dispersion  for  the  heavier  particles. 
These  figures  indicate  that  the  lighter  the  particles  the  greater  is 
their  horizontal  spread,  as  one  would  expect  In  the  case  a  B  5,  we 
have  also  included  results  for  t'  •»  10  and  t'  «=  10*°. 

2  A 

Taking  ~  400  cm  /sec  and  H  ~  20  km,  t'  =  10  corresponds  to 

6  26 
t  “  10  sec  ~  12  davs.  Assuming  n^  *  10  particles  and  z  «=  400  km, 

-4 

we  have  at  t'  »  10  the  following  table: 
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The  values  obtained  from  the  approximations  described  in  Section  2 
have  been  checked  at  several  points  by  direct  numerical  integration  of 
the  integral  N  using  an  IBM  1620  computer  for  a  -  1,  2,  5  and  10.  The 
agreement  of  the  two  has  been  found  to  be  satisfactory. 
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